Multibaryons with heavy flavors in the Skyrme model 
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We investigate the possible existence of multibaryons with heavy flavor quantum numbers using 
the bound state approach to the topological soliton model and the recently proposed approximation 
for multiskyrmion fields based on rational maps. We use an effective interaction lagrangian which 
consistently incorporates both chiral symmetry and the heavy quark symmetry including the cor- 

On i rections up to order O(l/mo). The model predicts some narrow heavy flavored multibaryon states 

ON ■ with baryon number four and seven. 

ON 

t-h ; 

PACS number(s): f2.39.Dc , 12.39. Fe , 14.20.Lq , 14.20.Mr 
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I. INTRODUCTION 

In recent years multibaryons have received a considerable amount of attention. A reason for this is that their 
hypothetical existence could provide valuable information about the nature of the strong interactions at low energies, 
t"**- ■ Perhaps the most celebrated example is that of the H-dibaryon predicted by Jaffe more than twenty years ago EJ. 

Since then the possible existence of some other exotic states has been investigated in various models. Of particular 
interest are those containing flavor (i.e. S , C , B ) quantum numbers. In fact, it has been speculated that strange 
matter could be stable Q . This has lead to numerous investigations of the properties of strange matter in bulk and 
. in finite lumps (for a recent review see Ref. ||). Moreover, with the advent of heavy ion colliders there is now the 
possibility of producing strange Q] and even charmed || multibaryonic states with rather low baryon number in the 
laboratory. These new developments provide further motivation for the study of the properties of multibaryons with 
heavy flavor quantum numbers. Most of the known predictions come from MIT bag model (see for example Refs. 
||) or non-relativistic quark model based calculations. Here, we will adopt a different point of view. We will asume 
^ that the heavy flavor multibaryons are formed by an SU{2) multiskyrmion with some heavy mesons bound to it. 
This is basically an extension of the bound state approach to strange hyperons originally introduced by Callan and 
Klebanov [Q] and later shown to describe heavier flavor baryons as well ||. A study of strange multibaryons within 
this approach has been recently presented in Ref. Q. There, a chiral lagrangian written in terms of pseudoscalar 
meson fields with some chiral symmetry breaking terms has been used as the effective lagrangian. As well known 
by now, although adequate for the light (up and down) and strange sectors, such type of effective lagr ang ian has to 
be modified when heavier flavors (e.g. charm) are incorporated. In that case, heavy quark symmetry pQ | has to be 
imposed. This symmetry requires that both the heavy pseudoscalar and the heavy vector fields appear explicitly in the 
effective lagrangian. Lagrangians which have both chiral symmetry and heavy quark symmetry have been described 
in the literature In our calculation we will adopt such type of lagrangian. As already mentioned, in our 

description the baryon number of the system comes from a non-trivial soliton configuration in the light sector. Until 
very recently only few multiskyrmion configurations (i.e. those with B < 4) were known. In 1997, however, after some 
demanding numerical work Battye and Sutcliff |lj| were able to identify those which are believed to be the lowest 
energy configurations with baryon number up to B = 9. Interestingly, all these configurations have the symmetries 
corresponding to the regular polyhedra. Even more important for our purposes, Houghton, Manton and Sutcliffe |P4| ] 
have exploited the similarities between the BPS monopoles and skyrmions to propose some ansatze based on rational 
maps. They have shown that such configurations approximate very well the numerically found lowest energy solutions 
with B < 9. In our investigations we will make use of these approximate ansatze. An interesting feature of these 
configurations is that for B > 1 the derivative of the radial soliton profile vanishes at the origin. Since to leading 
order in the inverse of the heavy quark mass uiq the heavy meson-soliton interaction is proportional to this quantity 
flEf we do not expect any bound state in that approximation. However, next-to- leading order corrections in 1/toq 
are required even to describe the spectrum of B = 1 heavy baryons |l6| ] . In the present work these corrections will be 
properly taken into account. 
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This article is organized as follows. In Sec. II we introduce the effective lagrangian together with the ansatz for the 
multiskyrmion configurations. In Sec. Ill we present our numerical results. Finally, in Sec. IV our main conclusions 
are given. 



II. THE LAGRANGIAN 



To describe the dynamics of the light and massive mesons interacting with each other we will consider an effective 
lagrangian of the following form p2| 



(1) 



where Ci is the effective light meson lagrangian. In the present work we will consider pions as the only explicit light 
degree of freedom and choose Ci to be simply the Skyrme lagrangian. Consequently, the effective lagrangian Ci , 
written in terms of the chiral field U — exp(ir • 7?// ff ), reads 



f r rr[d M rfd»U]+-^ i Tr 



Tr [U + C/ f - 2] 



(2) 



Here, f n is the pion decay constant and e is the so-called Skyrme parameter. In Eq.(l|), <j> represents the heavy 
pseudoscalar doublet and ip^ the corresponding vector doublet. For example, for charmed mesons 
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Moreover, / and g are the (pip^TT and ^^ip^n coupling constants, respectively, and 



Finally, in terms of the chiral field the currents and read 

a M = -(yWdpVU-VUdp</u?) . 



(3) 



(4) 
(5) 
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As usual in the bound state model we should first determine the static skyrmion background. For this purpose we 
introduce the rational map ansatz for the pion field. It reads 14 



Tr = f n F(r) n . (8) 
Here, F(r) is the (multi)skyrmion profile which depends on the radial coordinate only and n is a unit vector given by 

n = (2 t + 2 9f(i2) 3 + (1 - |i?| 2 ) k) , (9) 

where R is the rational map corresponding to a certain winding number B which is identified with the baryon number. 
Such map is usually written as a function of the complex variable z which is related to the usual spherical coordinates 
9, ip via stereographic projection, namely z = tan(0/2) exp(iip). For example, the map corresponding to the B = 1 
hedgehog ansatz is the identity map R — z. The explicit form of the maps corresponding to the other baryon numbers 
B < 9 can be found in Ref. p4| . Using Eq.(||) it is possible to obtain the expressions of the a M and currents to 
leading order in N c . The time components vanish at this order while the space components result 



= -- [F' T-hf J + sT-V l n] , 
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where we have introduced the short hand notation s = sinF, ^ = sin(-F/2). The radial profile function F(r) is 
determined by minimizing the soliton energy. Details of this procedure as well as plots of these functions for different 
baryon numbers are given in Ref. (lj] . 

To order N®, we have a system of heavy mesons moving in the static soliton background. To derive the explicit 
form of the relevant heavy meson-soliton lagrangian we need some consistent ansatze for the heavy meson fields. For 
the pseudoscalar field we use fl 



1 



47T 



(r, t) f ■ h x , 



(12) 



where ^ is a two-component spinor. To obtain the corresponding ansatz for the heavy vector meson field it is 
convenient to analize the coupling terms in the effective lagrangian. They are the last two terms in Eq.(^). From 
their structure and the form of the a M and currents in the static limit it is possible to see that the ansatz should 
have the form 



4-7T 
1 



W)\ (r, t) f % + i r ip2 (r, t) (n x V l n) • r] x 



(13) 
(14) 



Replacing Eqs.(|), (|T|), (|T|) and (jlj) in the effective la grangian, Eq.(|l|), we obtain that the heavy meson-soliton 
lagrangian L HM -soi is 



Lhm-soI — — 



dr r 2 I 



m 2 +2B^- 



+ ^fa + 2i^fa+fafa 



+2B 



+ 4^ - ( ^ + ^ - (- + 4- £<h) + K^-tbh 

r z \ r r J \ r r J r 



r 

-2gB 
+h.c. 



where s and £ have been already defined, (f = cos(F/2) and I is the angular integral 
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dVt (V l n a V l h a ) . 



(15) 



(16) 



The diagonalization of the hamiltonian obtained from Lhm-soI leads to a set of eigenvalue equations for the heavy 
meson field. They are 



u 4 + 2B^ T 



+ fB-fa = 
2 r 



(gF'f + up + g-) fa - 2gB-fa 2 = 
\ r I r 



(17) 
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ml-uj 2 + 2B^r fa + ~fF'<f> + — U ws - — ) 1 h ~ 2B —^2 = 



si 2 

2gB-^-fa - ufa 



+ + f^~<t> + - igF' + w) fa - g-lb'o -^-(F' + 2goj) fa - Ltf 
r 2r r r 2r r 



uj z - gujF' + 
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The numerical solution of this set of coupled equations supplemented with the appropiate boundary conditions 
provides the heavy meson energy ui for the different baryon numbers B. The corresponding results are discussed in 
the following section. 



III. NUMERICAL RESULTS 



In our numerical calculations we use two sets of parameters in the SU(2) sector. Set A corresponds to the case of 
massless pions and Set B to the case where the pion mass takes its empirical value — 138 MeV. In both cases, 
f v and e are adjusted so as to reproduce the empirical nucleon and A masses. The fitted values are f v — 64.5 MeV 
, e = 5.45 for Set A and f n — 54 MeV , e = 4.84 for Set B. With these parameters fixed and using for each baryon 
number B the rational map given in Ref. fb| ] , we obtain the profile F(r) that minimizes the mass of the soliton. 

We proceed to solve the bound state eigenvalue equations ([l7|)-(po|), using the values shown in Table I for the 
parameters that appear in the heavy meson lagrangian, Eq. ([!]). For the pseudoscalar and vector meson masses we 
use the empirical values. On the other hand, since little is known about the heavy meson coupling constants, we use 
the heavy quark symmetry relation jl2| 

/ = 2 g (21) 

as a guideline in order to estimate /. We take g to be the value given by the non-relativistic quark model: g = —0.75. 
As discussed in a recent analysis fl8l , this value is compatible with the upper bound g 2 < .5 established by the 
experimental upper limit for the decay width T(D* + ) < 131 keV set by the ACCMOR Collaboration Q. 
Our results for the heavy meson binding energies Eb defined by 

£b = - lo b (22) 

are shown in Table II. Also listed in Table II are the soliton masses per baryon number M so i taken from Ref. 0. 

The general structure shown by the binding energies is in qualitative agreement with what was found in Ref. |{J 
for the strangeness case. For the heavy flavors we also find that the binding energy decreases with increasing baryon 
number, except for the crossings that take place at B = 4 and B = 7. This general behaviour is the opposite to what 
was found in Ref. p7| . It should be pointed out that, due to the absence of explicit vector mesons in the corresponding 
effective action, the approach used in |lj] is expected to be less accurate than the one followed in the present work. 

From Table II we also notice that, for heavy flavors, the binding is stronger for Set A than for Set B. This can 
be understood as follows. In the wiq — > oo limit the heavy flavored meson would be concentrated at the origin [ jl5| , 
wrapped by the soliton. Thus, it would only probe the potential at this point. For B = 1 such potential is basically 
proportional to ^^'(O)! and attractive. As well known (see, e.g., Fig.l of Ref. ^0|), |-F"(0)| is larger in the massless 
case. That leads to the observed behaviour. Similar analysis can be done for higher values of B. For the strangeness 
case the behaviour of the profile function at medium distances becomes important and the situation is reversed || . 

In order to study the stability of the heavy multiskyrmions we will only consider the mass of the background 
multiskyrmion and the binding energy of the bound mesons. This should be a good approximation since the non- 
adiabatic corrections are expected to be small as a consequence of the rather large values of the moments of inertia 
involved pl[ . In the following we will focus on those states which have flavor number equal to their baryon number. 
These states are of particular interest since in the strange sector the analog states, namely those with Y = 0, have 
been found to be stable for some values of B (2^|J . 

Using 

I B = M 1 +M B -i-M B (23) 

with Mb = B(M so i + m^, — eb) and the values given in Table II, we get the ionization energies Ib shown in Table 
III. We observe that the only heavy flavored states that may be stable are those with B = 4 and B = 7. In Table 
IV we summarize the energies for the other possible strong decays of these states. We observe that, although these 
states are stable against strong decays into two fragments, some decays into a larger number of fragments are allowed. 
However, since the usual phase space factors tend to suppress the decay rates as the number of fragments n in the 
final state increases we expect them to be quite narrow. For instance, in the case of the charmed heptalambda the 
decay width will be very small since the only allowed decay mode is the one that has seven A c in the final state (Set 
A). 

Therefore the present model predicts, both for charm and bottom, narrow heavy multibaryon states with baryon 
number four and seven. The main reason for the unstability of these particles can be traced back to the rather large 
difference that exists between the B — 1 and B > 1 binding energies. This can be easily understood noting that, 
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to leading order in l/mg, the only non- vanishing binding energy would be that corresponding to B = 1 since for 
B > 1 the radial derivative of the soliton profile function vanishes at the origin. On the other hand, in the case of 
strangeness the meson wavefunction is wider and, therefore, much less sensitive to the value of the potential at the 
origin. Consequently, the gap between the binding energies of B = 1 and B > 1 is much smaller and the corresponding 
multibaryon states with baryon number four and seven turn out to be absolutely stable against strong decays M . 



IV. CONCLUSIONS 



In this work we have studied the masses of heavy multibaryon configurations using the bound state approach to 
the Skyrme model. This is a natural extension of previous work done in the strange sector. In order to consider the 
heavier flavors of charm and bottom in a consistent way, however, it is important to take into account the heavy quark 
symmetry. This is accomplished by the lagrangian given in Ref . |l2| . An important feature of this effective heavy 
meson lagrangian is that it contains the degrees of freedom of the heavy scalar meson and the heavy vector meson 
explicitly, which leads to a system of four coupled equations for the bound state problem instead of just one equation 
as in Ref. The baryon number is carried by the soliton configuration of the light background fields, for which we 
have used the expressions in terms of the rational maps given in Ref. Jl4|] . 

We obtained solutions for the bound state equations for B < 9. As in the strangeness case we find that the binding 
energy decreases with the baryon number and that the B = 4 and B — 7 states are the most stable against strong 
decays. However, for the charm and bottom flavors these states are not absolutely stable. Still, they are expected 
to be quite narrow since only decays into final states with three or more fragments are energetically allowed. We do 
not expect that the collective quantization of the soliton-meson bound system will change this picture. It would be 
important, however, to estimate the zero point energy |p3| of these multiskyrmion configurations. As discussed in 
Ref. J24| this contribution may be the cause for the H particle, which appears almost at threshold, to be unbound. 
In any case, since the predicted tetralambda and heptalambda are more strongly bound against two particle decays 
than the H, this contribution is not expected to be so important so as to open those leading decay channels. 
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TABLE I. Masses and coupling constants for the heavy meson lagrangian. 







charm 




bottom 






1867 MeV 




5279 MeV 






2010 MeV 




5325 MeV 


f 




-3016 MeV 




-7988 MeV 


g 




-0.75 




-0.75 


TABLE II. 


Meson binding energies eb (in MeV) for the case of massless pions (Set A) and 


massive pions (Set B) as a 


function of the baryon number B. 


Also listed are the soliton masses per baryon unit M ao i (in MeV) which are taken from Ref. 


[9]. 












Msol 


e b (charm) 




£b (bottom) 


B 


Set A 


Set B Set A Set B 




Set A Set B 


1 


863 


864 383 328 




554 474 


2 


847 


848 321 272 




438 374 


3 


830 


832 300 255 




408 351 


4 


797 


798 301 256 




407 350 


5 


804 


808 287 245 




391 339 


6 


797 


802 283 243 




386 336 


7 


776 


780 288 247 




392 341 


8 


784 


790 280 242 




383 335 


9 


787 


796 275 239 




378 332 


TABLE III. 


Ionization energies 


1b (in MeV) of the charm and bottom multilambdas 


in the case of massless pions (Set A) 


and massive pions (Set B). 












Ib (charm ) 




Is (bottom ) 


B 


Set A 


Set B 


Set A 


Set B 


2 


-92 


-80 


-200 


-168 


3 


-58 


-43 


-139 


-105 


4 


86 


99 


15 


41 


5 


-121 


-111 


-196 


-163 


6 


-19 


-3 


-92 


-61 


7 


148 


159 


87 


113 


8 


-136 


-117 


-211 


-177 


9 


-96 


-93 


-164 


-146 
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TABLE IV. Energy balance (in MeV) for the multilambda strong decays in the case of massless pions (Set A) and massive 
pions (Set B). The first column indicates the number of fragments n in the final state. 









f V| Q T"TTI 
I^IICLI 111 




UU Li LLIIll 










Set A 


Set B 


Set A 


Set B 


2 


M 4 A 


- 2M 2 a 


-120 


-136 


-76 


-104 


3 


M 4 a 


- (2M A + M 2 a) 


-28 


-56 


124 


64 


4 


M 4 A 


- 4Ma 


64 


24 


324 


232 


2 


M 7A 


- (M 2A + Ms a) 


-221 


-236 


-195 


-220 




M 7A 


- (M 3A + M 4A ) 


-158 


-168 


-138 


-162 


3 


Mr A 
ivi 7 a 


— f Ma 4- Mo a 4- M/i a 1 


-100 


-125 


1 


-57 




M 7A 


- (M A + 2M 3A ) 


-244 


-267 


-153 


-203 




M 7 A 


- (2M A + M 5A ) 


-129 


-156 


5 


-52 




M 7A 


- (2M 2A + M 3A ) 


-278 


-304 


-214 


-266 


4 


M 7A 


- (M A + 3M 2A ) 


-220 


-261 


-75 


-161 




M 7A 


- (2M A + M 2A + M 3A ) 


-186 


-224 


-14 


-98 




M 7A 


- (3M A + M 4A ) 


-8 


-45 


201 


111 


5 


M 7A 


-(3M A + 2M 2A ) 


-128 


-181 


125 


7 




M 7A 


- (4M A + M 3 a) 


-94 


-144 


186 


70 


6 


M 7A 


- (5M A + M 2A ) 


-36 


-101 


325 


175 


7 


M 7A 


- (7M A ) 


56 


-21 


525 


343 
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